Incomplete lineage sorting and phenotypic evolution in marsupials.

Incomplete lineage sorting (ILS) makes ancestral genetic polymorphisms persist during rapid speciation events, inducing incongruences between gene trees and species trees. ILS has complicated phylogenetic inference in many lineages, including hominids. However, we lack empirical evidence that ILS leads to incongruent phenotypic variation. Here, we performed phylogenomic analyses to show that the South American monito del monte is the sister lineage of all Australian marsupials, although over 31% of its genome is closer to the Diprotodontia than to other Australian groups due to ILS during ancient radiation. Pervasive conflicting phylogenetic signals across the whole genome are consistent with some of the morphological variation among extant marsupials. We detected hundreds of genes that experienced stochastic fixation during ILS, encoding the same amino acids in non-sister species. Using functional experiments, we confirm how ILS may have directly contributed to hemiplasy in morphological traits that were established during rapid marsupial speciation ca. 60 mya.

Tracing the Enterococci from Paleozoic Origins to the Hospital.

We examined the evolutionary history of leading multidrug resistant hospital pathogens, the enterococci, to their origin hundreds of millions of years ago. Our goal was to understand why, among the vast diversity of gut flora, enterococci are so well adapted to the modern hospital environment. Molecular clock estimation, together with analysis of their environmental distribution, phenotypic diversity, and concordance with host fossil records, place the origins of the enterococci around the time of animal terrestrialization, 425-500 mya. Speciation appears to parallel the diversification of hosts, including the rapid emergence of new enterococcal species following the End Permian Extinction. Major drivers of speciation include changing carbohydrate availability in the host gut. Life on land would have selected for the precise traits that now allow pathogenic enterococci to survive desiccation, starvation, and disinfection in the modern hospital, foreordaining their emergence as leading hospital pathogens.

Mitochondrial DNA variation in human radiation and disease.

Environmental adaptation, predisposition to common diseases, and, potentially, speciation may all be linked through the adaptive potential of mitochondrial DNA (mtDNA) alterations of bioenergetics. This Perspective synthesizes evidence that human mtDNA variants may be adaptive or deleterious depending on environmental context and proposes that the accrual of mtDNA variation could contribute to animal speciation via adaptation to marginal environments.

Cenozoic history of the tropical marine biodiversity hotspot.

The region with the highest marine biodiversity on our planet is known as the Coral Triangle or Indo-Australian Archipelago (IAA)1,2. Its enormous biodiversity has long attracted the interest of biologists; however, the detailed evolutionary history of the IAA biodiversity hotspot remains poorly understood3. Here we present a high-resolution reconstruction of the Cenozoic diversity history of the IAA by inferring speciation-extinction dynamics using a comprehensive fossil dataset. We found that the IAA has exhibited a unidirectional diversification trend since about 25 million years ago, following a roughly logistic increase until a diversity plateau beginning about 2.6 million years ago. The growth of diversity was primarily controlled by diversity dependency and habitat size, and also facilitated by the alleviation of thermal stress after 13.9 million years ago. Distinct net diversification peaks were recorded at about 25, 20, 16, 12 and 5 million years ago, which were probably related to major tectonic events in addition to climate transitions. Key biogeographic processes had far-reaching effects on the IAA diversity as shown by the long-term waning of the Tethyan descendants versus the waxing of cosmopolitan and IAA taxa. Finally, it seems that the absence of major extinctions and the Cenozoic cooling have been essential in making the IAA the richest marine biodiversity hotspot on Earth.

A lethal mitonuclear incompatibility in complex I of natural hybrids.

The evolution of reproductive barriers is the first step in the formation of new species and can help us understand the diversification of life on Earth. These reproductive barriers often take the form of hybrid incompatibilities, in which alleles derived from two different species no longer interact properly in hybrids1-3. Theory predicts that hybrid incompatibilities may be more likely to arise at rapidly evolving genes4-6 and that incompatibilities involving multiple genes should be common7,8, but there has been sparse empirical data to evaluate these predictions. Here we describe a mitonuclear incompatibility involving three genes whose protein products are in physical contact within respiratory complex I of naturally hybridizing swordtail fish species. Individuals homozygous for mismatched protein combinations do not complete embryonic development or die as juveniles, whereas those heterozygous for the incompatibility have reduced complex I function and unbalanced representation of parental alleles in the mitochondrial proteome. We find that the effects of different genetic interactions on survival are non-additive, highlighting subtle complexity in the genetic architecture of hybrid incompatibilities. Finally, we document the evolutionary history of the genes involved, showing signals of accelerated evolution and evidence that an incompatibility has been transferred between species via hybridization.

Hybrid speciation driven by multilocus introgression of ecological traits.

Hybridization allows adaptations to be shared among lineages and may trigger the evolution of new species1,2. However, convincing examples of homoploid hybrid speciation remain rare because it is challenging to demonstrate that hybridization was crucial in generating reproductive isolation3. Here we combine population genomic analysis with quantitative trait locus mapping of species-specific traits to examine a case of hybrid speciation in Heliconius butterflies. We show that Heliconius elevatus is a hybrid species that is sympatric with both parents and has persisted as an independently evolving lineage for at least 180,000 years. This is despite pervasive and ongoing gene flow with one parent, Heliconius pardalinus, which homogenizes 99% of their genomes. The remaining 1% introgressed from the other parent, Heliconius melpomene, and is scattered widely across the H. elevatus genome in islands of divergence from H. pardalinus. These islands contain multiple traits that are under disruptive selection, including colour pattern, wing shape, host plant preference, sex pheromones and mate choice. Collectively, these traits place H. elevatus on its own adaptive peak and permit coexistence with both parents. Our results show that speciation was driven by introgression of ecological traits, and that speciation with gene flow is possible with a multilocus genetic architecture.

A continuous fish fossil record reveals key insights into adaptive radiation.

Adaptive radiations have been instrumental in generating a considerable amount of life's diversity. Ecological opportunity is thought to be a prerequisite for adaptive radiation1, but little is known about the relative importance of species' ecological versatility versus effects of arrival order in determining which lineage radiates2. Palaeontological records that could help answer this are scarce. In Lake Victoria, a large adaptive radiation of cichlid fishes evolved in an exceptionally short and recent time interval3. We present a rich continuous fossil record extracted from a series of long sediment cores along an onshore-offshore gradient. We reconstruct the temporal sequence of events in the assembly of the fish community from thousands of tooth fossils. We reveal arrival order, relative abundance and habitat occupation of all major fish lineages in the system. We show that all major taxa arrived simultaneously as soon as the modern lake began to form. There is no evidence of the radiating haplochromine cichlid lineage arriving before others, nor of their numerical dominance upon colonization; therefore, there is no support for ecological priority effects. However, although many taxa colonized the lake early and several became abundant, only cichlids persisted in the new deep and open-water habitats once these emerged. Because these habitat gradients are also known to have played a major role in speciation, our findings are consistent with the hypothesis that ecological versatility was key to adaptive radiation, not priority by arrival order nor initial numerical dominance.

Assembly of functional diversity in an oceanic island flora.

Oceanic island floras are well known for their morphological peculiarities and exhibit striking examples of trait evolution1-3. These morphological shifts are commonly attributed to insularity and are thought to be shaped by the biogeographical processes and evolutionary histories of oceanic islands2,4. However, the mechanisms through which biogeography and evolution have shaped the distribution and diversity of plant functional traits remain unclear5. Here we describe the functional trait space of the native flora of an oceanic island (Tenerife, Canary Islands, Spain) using extensive field and laboratory measurements, and relate it to global trade-offs in ecological strategies. We find that the island trait space exhibits a remarkable functional richness but that most plants are concentrated around a functional hotspot dominated by shrubs with a conservative life-history strategy. By dividing the island flora into species groups associated with distinct biogeographical distributions and diversification histories, our results also suggest that colonization via long-distance dispersal and the interplay between inter-island dispersal and archipelago-level speciation processes drive functional divergence and trait space expansion. Contrary to our expectations, speciation via cladogenesis has led to functional convergence, and therefore only contributes marginally to functional diversity by densely packing trait space around shrubs. By combining biogeography, ecology and evolution, our approach opens new avenues for trait-based insights into how dispersal, speciation and persistence shape the assembly of entire native island floras.

Origination of the modern-style diversity gradient 15 million years ago.

The latitudinal diversity gradient (LDG) is a prevalent feature of modern ecosystems across diverse clades1-4. Recognized for well over a century, the causal mechanisms for LDGs remain disputed, in part because numerous putative drivers simultaneously covary with latitude1,3,5. The past provides the opportunity to disentangle LDG mechanisms because the relationships among biodiversity, latitude and possible causal factors have varied over time6-9. Here we quantify the emergence of the LDG in planktonic foraminifera at high spatiotemporal resolution over the past 40 million years, finding that a modern-style gradient arose only 15 million years ago. Spatial and temporal models suggest that LDGs for planktonic foraminifera may be controlled by the physical structure of the water column. Steepening of the latitudinal temperature gradient over 15 million years ago, associated with an increased vertical temperature gradient at low latitudes, may have enhanced niche partitioning and provided more opportunities for speciation at low latitudes. Supporting this hypothesis, we find that higher rates of low-latitude speciation steepened the diversity gradient, consistent with spatiotemporal patterns of depth partitioning by planktonic foraminifera. Extirpation of species from high latitudes also strengthened the LDG, but this effect tended to be weaker than speciation. Our results provide a step change in understanding the evolution of marine LDGs over long timescales.

Genetic and neural mechanisms that inhibit Drosophila from mating with other species.

Genetically hard-wired neural mechanisms must enforce behavioral reproductive isolation because interspecies courtship is rare even in sexually naïve animals of most species. We find that the chemoreceptor Gr32a inhibits male D. melanogaster from courting diverse fruit fly species. Gr32a recognizes nonvolatile aversive cues present on these reproductively dead-end targets, and activity of Gr32a neurons is necessary and sufficient to inhibit interspecies courtship. Male-specific Fruitless (Fru(M)), a master regulator of courtship, also inhibits interspecies courtship. Gr32a and Fru(M) are not coexpressed, but Fru(M) neurons contact Gr32a neurons, suggesting that these genes influence a shared neural circuit that inhibits interspecies courtship. Gr32a and Fru(M) also suppress within-species intermale courtship, but we show that distinct mechanisms preclude sexual displays toward conspecific males and other species. Although this chemosensory pathway does not inhibit interspecies mating in D. melanogaster females, similar mechanisms appear to inhibit this behavior in many other male drosophilids.

Diverse pathways to speciation revealed by marine snails.

Speciation is a key evolutionary process that is not yet fully understood. Combining population genomic and ecological data from multiple diverging pairs of marine snails (Littorina) supports the search for speciation mechanisms. Placing pairs on a one-dimensional speciation continuum, from undifferentiated populations to species, obscured the complexity of speciation. Adding multiple axes helped to describe either speciation routes or reproductive isolation in the snails. Divergent ecological selection repeatedly generated barriers between ecotypes, but appeared less important in completing speciation while genetic incompatibilities played a key role. Chromosomal inversions contributed to genomic barriers, but with variable impact. A multidimensional (hypercube) approach supported framing of questions and identification of knowledge gaps and can be useful to understand speciation in many other systems.

A three-dimensional genomics view for speciation research.

Genomic information is folded in a three-dimensional (3D) structure, a rarely explored evolutionary driver of speciation. Technological advances now enable the study of 3D genome structures (3DGSs) across the Tree of Life. At the onset of 3D speciation genomics, we discuss the putative roles of 3DGSs in speciation.

Drivers and dynamics of a massive adaptive radiation in cichlid fishes.

Adaptive radiation is the likely source of much of the ecological and morphological diversity of life1-4. How adaptive radiations proceed and what determines their extent remains unclear in most cases1,4. Here we report the in-depth examination of the spectacular adaptive radiation of cichlid fishes in Lake Tanganyika. On the basis of whole-genome phylogenetic analyses, multivariate morphological measurements of three ecologically relevant trait complexes (body shape, upper oral jaw morphology and lower pharyngeal jaw shape), scoring of pigmentation patterns and approximations of the ecology of nearly all of the approximately 240 cichlid species endemic to Lake Tanganyika, we show that the radiation occurred within the confines of the lake and that morphological diversification proceeded in consecutive trait-specific pulses of rapid morphospace expansion. We provide empirical support for two theoretical predictions of how adaptive radiations proceed, the 'early-burst' scenario1,5 (for body shape) and the stages model1,6,7 (for all traits investigated). Through the analysis of two genomes per species and by taking advantage of the uneven distribution of species in subclades of the radiation, we further show that species richness scales positively with per-individual heterozygosity, but is not correlated with transposable element content, number of gene duplications or genome-wide levels of selection in coding sequences.

Genetic drift promotes and recombination hinders speciation on holey fitness landscapes.

Dobzhansky and Muller proposed a general mechanism through which microevolution, the substitution of alleles within populations, can cause the evolution of reproductive isolation between populations and, therefore, macroevolution. As allopatric populations diverge, many combinations of alleles differing between them have not been tested by natural selection and may thus be incompatible. Such genetic incompatibilities often cause low fitness in hybrids between species. Furthermore, the number of incompatibilities grows with the genetic distance between diverging populations. However, what determines the rate and pattern of accumulation of incompatibilities remains unclear. We investigate this question by simulating evolution on holey fitness landscapes on which genetic incompatibilities can be identified unambiguously. We find that genetic incompatibilities accumulate more slowly among genetically robust populations and identify two determinants of the accumulation rate: recombination rate and population size. In large populations with abundant genetic variation, recombination selects for increased genetic robustness and, consequently, incompatibilities accumulate more slowly. In small populations, genetic drift interferes with this process and promotes the accumulation of genetic incompatibilities. Our results suggest a novel mechanism by which genetic drift promotes and recombination hinders speciation.

Chromosomal evolution, environmental heterogeneity, and migration drive spatial patterns of species richness in Calochortus (Liliaceae).

We used nuclear genomic data and statistical models to evaluate the ecological and evolutionary processes shaping spatial variation in species richness in Calochortus (Liliaceae, 74 spp.). Calochortus occupies diverse habitats in the western United States and Mexico and has a center of diversity in the California Floristic Province, marked by multiple orogenies, winter rainfall, and highly divergent climates and substrates (including serpentine). We used sequences of 294 low-copy nuclear loci to produce a time-calibrated phylogeny, estimate historical biogeography, and test hypotheses regarding drivers of present-day spatial patterns in species number. Speciation and species coexistence require reproductive isolation and ecological divergence, so we examined the roles of chromosome number, environmental heterogeneity, and migration in shaping local species richness. Six major clades-inhabiting different geographic/climatic areas, and often marked by different base chromosome numbers (n = 6 to 10)-began diverging from each other ~10.3 Mya. As predicted, local species number increased significantly with local heterogeneity in chromosome number, elevation, soil characteristics, and serpentine presence. Species richness is greatest in the Transverse/Peninsular Ranges where clades with different chromosome numbers overlap, topographic complexity provides diverse conditions over short distances, and several physiographic provinces meet allowing immigration by several clades. Recently diverged sister-species pairs generally have peri-patric distributions, and maximum geographic overlap between species increases over the first million years since divergence, suggesting that chromosomal evolution, genetic divergence leading to gametic isolation or hybrid inviability/sterility, and/or ecological divergence over small spatial scales may permit species co-occurrence.

Rapid vertebrate speciation via isolation, bottlenecks, and drift.

Speciation is often driven by selective processes like those associated with viability, mate choice, or local adaptation, and "speciation genes" have been identified in many eukaryotic lineages. In contrast, neutral processes are rarely considered as the primary drivers of speciation, especially over short evolutionary timeframes. Here, we describe a rapid vertebrate speciation event driven primarily by genetic drift. The White Sands pupfish (Cyprinodon tularosa) is endemic to New Mexico's Tularosa Basin where the species is currently managed as two Evolutionarily significant units (ESUs) and is of international conservation concern (Endangered). Whole-genome resequencing data from each ESU showed remarkably high and uniform levels of differentiation across the entire genome (global FST ≈ 0.40). Despite inhabiting ecologically dissimilar springs and streams, our whole-genome analysis revealed no discrete islands of divergence indicative of strong selection, even when we focused on an array of candidate genes. Demographic modeling of the joint allele frequency spectrum indicates the two ESUs split only ~4 to 5 kya and that both ESUs have undergone major bottlenecks within the last 2.5 millennia. Our results indicate the genome-wide disparities between the two ESUs are not driven by divergent selection but by neutral drift due to small population sizes, geographic isolation, and repeated bottlenecks. While rapid speciation is often driven by natural or sexual selection, here we show that isolation and drift have led to speciation within a few thousand generations. We discuss these evolutionary insights in light of the conservation management challenges they pose.

Extreme elevational migration spurred cryptic speciation in giant hummingbirds.

The ecoevolutionary drivers of species niche expansion or contraction are critical for biodiversity but challenging to infer. Niche expansion may be promoted by local adaptation or constrained by physiological performance trade-offs. For birds, evolutionary shifts in migratory behavior permit the broadening of the climatic niche by expansion into varied, seasonal environments. Broader niches can be short-lived if diversifying selection and geography promote speciation and niche subdivision across climatic gradients. To illuminate niche breadth dynamics, we can ask how "outlier" species defy constraints. Of the 363 hummingbird species, the giant hummingbird (Patagona gigas) has the broadest climatic niche by a large margin. To test the roles of migratory behavior, performance trade-offs, and genetic structure in maintaining its exceptional niche breadth, we studied its movements, respiratory traits, and population genomics. Satellite and light-level geolocator tracks revealed an >8,300-km loop migration over the Central Andean Plateau. This migration included a 3-wk, ~4,100-m ascent punctuated by upward bursts and pauses, resembling the acclimatization routines of human mountain climbers, and accompanied by surging blood-hemoglobin concentrations. Extreme migration was accompanied by deep genomic divergence from high-elevation resident populations, with decisive postzygotic barriers to gene flow. The two forms occur side-by-side but differ almost imperceptibly in size, plumage, and respiratory traits. The high-elevation resident taxon is the world's largest hummingbird, a previously undiscovered species that we describe and name here. The giant hummingbirds demonstrate evolutionary limits on niche breadth: when the ancestral niche expanded due to evolution (or loss) of an extreme migratory behavior, speciation followed.

TRAILS: Tree reconstruction of ancestry using incomplete lineage sorting.

Genome-wide genealogies of multiple species carry detailed information about demographic and selection processes on individual branches of the phylogeny. Here, we introduce TRAILS, a hidden Markov model that accurately infers time-resolved population genetics parameters, such as ancestral effective population sizes and speciation times, for ancestral branches using a multi-species alignment of three species and an outgroup. TRAILS leverages the information contained in incomplete lineage sorting fragments by modelling genealogies along the genome as rooted three-leaved trees, each with a topology and two coalescent events happening in discretized time intervals within the phylogeny. Posterior decoding of the hidden Markov model can be used to infer the ancestral recombination graph for the alignment and details on demographic changes within a branch. Since TRAILS performs posterior decoding at the base-pair level, genome-wide scans based on the posterior probabilities can be devised to detect deviations from neutrality. Using TRAILS on a human-chimp-gorilla-orangutan alignment, we recover speciation parameters and extract information about the topology and coalescent times at high resolution.

Recent secondary contact, genome-wide admixture, and asymmetric introgression of neo-sex chromosomes between two Pacific island bird species.

Secondary contact between closely related taxa represents a "moment of truth" for speciation-an opportunity to test the efficacy of reproductive isolation that evolved in allopatry and to identify the genetic, behavioral, and/or ecological barriers that separate species in sympatry. Sex chromosomes are known to rapidly accumulate differences between species, an effect that may be exacerbated for neo-sex chromosomes that are transitioning from autosomal to sex-specific inheritance. Here we report that, in the Solomon Islands, two closely related bird species in the honeyeater family-Myzomela cardinalis and Myzomela tristrami-carry neo-sex chromosomes and have come into recent secondary contact after ~1.1 my of geographic isolation. Hybrids of the two species were first observed in sympatry ~100 years ago. To determine the genetic consequences of hybridization, we use population genomic analyses of individuals sampled in allopatry and in sympatry to characterize gene flow in the contact zone. Using genome-wide estimates of diversity, differentiation, and divergence, we find that the degree and direction of introgression varies dramatically across the genome. For sympatric birds, autosomal introgression is bidirectional, with phenotypic hybrids and phenotypic parentals of both species showing admixed ancestry. In other regions of the genome, however, the story is different. While introgression on the Z/neo-Z-linked sequence is limited, introgression of W/neo-W regions and mitochondrial sequence (mtDNA) is highly asymmetric, moving only from the invading M. cardinalis to the resident M. tristrami. The recent hybridization between these species has thus enabled gene flow in some genomic regions but the interaction of admixture, asymmetric mate choice, and/or natural selection has led to the variation in the amount and direction of gene flow at sex-linked regions of the genome.

Mapping seasonal migration in a songbird hybrid zone -- heritability, genetic correlations, and genomic patterns linked to speciation.

Seasonal migration is a widespread behavior relevant for adaptation and speciation, yet knowledge of its genetic basis is limited. We leveraged advances in tracking and sequencing technologies to bridge this gap in a well-characterized hybrid zone between songbirds that differ in migratory behavior. Migration requires the coordinated action of many traits, including orientation, timing, and wing morphology. We used genetic mapping to show these traits are highly heritable and genetically correlated, explaining how migration has evolved so rapidly in the past and suggesting future responses to climate change may be possible. Many of these traits mapped to the same genomic regions and small structural variants indicating the same, or tightly linked, genes underlie them. Analyses integrating transcriptomic data indicate cholinergic receptors could control multiple traits. Furthermore, analyses integrating genomic differentiation further suggested genes underlying migratory traits help maintain reproductive isolation in this hybrid zone.